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TECHNICAL MEMORANDUM 1264 


CALCULATION OF TEE LATERAL-DYNAMIC STABILITY OF AIRCRAFT* 

By A. Raikh 

The object of this report is to present a method of computing 
the lateral-dynamic stability of an airplane. 

Graphs and formulas are given with the aid of which all the 
aerodynamic coefficients required for computing the lateral dynamic 
stability can be determined. A number of numerical examples are 
given for obtaining the stability derivatives and solving the 
characteristic-stability equation. Approximate formulas are derived 
with the aid of which rapid preliminary computations may be made and 
the stability coefficients corrected for certain modifications of 
the airplane. A derivation of the lateral-dynamic-stability eq.ua- 
tions is included. 


INTRODUCTION h ■ 

In the present .stage of development of airplane design, con- 
stantly increasing requirements are imposed on the airplane as 
regards stability and maneuverability. If the area of the vertical 
tail surface and the transverse dihedral of the wing are unfavorably 
chosen, the airplane will be subjected to large lateral motions at 
the least gust of wind and the pilot will be under the constant 
necessity of applying the controls. Some airplanes because of their 
flying qualities are generally approved by pilots. Computations show 
that such airplanes are spirally unstable only at sufficiently large 
angles of attack but are spirally stable at small angles of attack. 

A method of computation of the lateral- dynamic stability is 
presented herein. The computation is illustrated by examples of 
the computation of a Northrup 2E airplane, which has received the 
highest commendation of pilots for its piloting characteristics with 
regard to stability and maneuverability. 

The difficulty of computing the dynamic stability lies in the 
fact that a large number of aerodynamic coefficients, the so-called 
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rotary derivatives, are required for the computation. These deriva- 
tives are not encountered in the usual aerodynamic computation. In 
order to facilitate the computation, a number of curves are given 
with the aid of which these magnitudes are rapidly and simply deter- 
mined. The numerical examples given for computing the coefficients 
further simplify the problem. 

A comparison of the aerodynamio coefficient obtained by compu- 
tation from graphs and that determined from tests in the wind tunnel 
is presented herein. In order to oheck the accuracy of the final 
data, a comparison is given of the period of oscillation and the 
time taken for the damping of the lateral motion as obtained from 
flight tests and by computation, respectively. The comparison shows 
that the computation gives good results. 

Notwithstanding the apparent complexity, the entire computation 
on the lateral stability from the available aerodynamic coefficients 
takes no longer than 1.5 to 2 hours. The author expresses his thanks 
to A. I. Silman, a candidate for a technical degree, for a number of 
valuable suggestions utilised in this work. 


SYSTEM OF AXES AND NOTATION 


The system of axes is shown in figure 1. The origin of the 
coordinates is at the center of gravity of the airplane. 

The X-axis is in the plane of symmetry of the airplane and dir- 
ected parallel to the velocity in steady flight. The Y-axis is per- 
pendicular to the plane of symmetry of the airplane and directed to 
the right of the pilot. The Z-axis is in the plane of symmetry per- 
pendicular to the X-axis and directed downwards. 

X, Y, Z forces along corresponding axes; positive if 

direction agrees with direction of axes 


L, M, N 


u, v, w 


P, r 


moments' about corresponding axes: moment L is 
positive if it causes right wing to drop; 
moment M is positive if it causes tail to 
drop; moment N is positive if it causes 
backward motion of right wing 

projections of linear velocity on corresponding 
axes; positive if direction of projections 
agrees with directions of axes ^ 

projections of angular velocities; have same posi- 
tive directions as moments L, M, N 
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A, B, C 
E 
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X 


x 


S 



b 


angles mad© "by given direction with. X-, Y-, Z- 
axes respectively; have same positive direc- 
tions as L, M, N 

moments of inertia about X-, Y-, Z-axes 

the centrifugal moment of inertia [Ed. note, 
product of inertia with respect to X- and Y- 
axis] 

angle of attack, degrees 

angle of attack computed from line of zero lift, 
degrees 

sideslip angle, £3 = -arc sin v/V, radians (jSACA 
comment: {3 as used in this report is of 

opposite sign to p as used in American and 
British reports]) 

angle of inclination of flight path to horizon- 
tal (positive for lift), degrees 

angle of dihedral of wing, angle between plane 
of chords and plane at right angles to plane 
of symmetry and passing through chord at tip, 
degrees 

angle of sveepback, angle between focal line at 
0.25 chord computed from leading edge and 
plane perpendicular to axis of fuselage, 
degrees 

angle between line joining center of gravity 
of airplane with geometric center of verti- 
cal tail surface (with part of fuselage) 
and line of zero lift, degrees 

area of wing 

area of vertical tail with part of fuselage 
(fig. 11) 

area of lateral projection of fuselage 


span 
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c 0 
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A . b 2 /S 

T •= °o/ c f 

l nf 

*f 

P 

*vt 





y P 


d°y 

^3 


. ^ 
Z P B W 


chord at root of wing 

chord at tip of ving, effective tip chord is 
determined by prolonging lines of leading 
and. trailing edges 

aspect ratio 

taper of wing 

distance from nose of fuselage to center of 
gravity of airplane 

distance from rudder hinge to center of gravity 
length of fuselage 
density of the air 

coefficient of effectiveness of vertical tail 
surface 

coefficient of lateral force 


coefficient of rolling moment 


coefficient of yawing moment 


lift force coefficient 


nondimensional derivative of lateral force with 
respect to sideslip angle (3 in radians) 

nondimensional derivative of rolling moment with 
respect to sideslip angle 
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nondimensional derivative of yawing moment with, 
respect to sideslip angle 

nondimensional derivative of rolling moment with, 
respect to angular' velocity of roll 


nondimensional derivative of yawing moment with 
respect to angular velocity of roll 


nondimensional derivative of rolling moment with 
respect to angular velocity of yaw 


nondimensional derivative of yawing moment with 
respect to angular velocity of yaw 


coefficient of moment of inertia about X-axis 


coefficient of moment of inertia about Z-axis 

weight of airplane 
mass of airplane 


relative density of airplane 


T * 


m 

pSV 


unit of time 


t 



nondimensional time 


a 4> a 3 } Ag, A Aq coefficients of characteristic equation 
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1. EQUATIONS OF LATERAL-DYNAMIC STABILITY OF AIRPLANE 

If the rectilinear steady flight of an airplane is considered, 
the motion of airplane after a certain disturbance is described by 
a system of sir equations with six variables: 

u, w, q, and p, p, r 

If. only small deviations of the airplane from steady recti- 
linear flight are considered, it may be shown (reference l) t ha t 
the system of six differential equations breaks up into two inde- 
pendent systems of equations. In the first system of differential 
equations the following variables enter: 

u projection of linear velocity on X-axis 

w projection of linear velocity on Z-axis 

q projection of angular velooity on Y-axis 

In the first system of equations, there enter only variables 
characterizing the motion of the airplane in the plane of symmetry 
(that is in the plane XOZ) . In the second group of equations there 
enter the variables: 

P angle of sideslip, determined by projection of velocity on 
Y-axis (sideslip velocity v) 

p projection of angular velocity on X-axis 

r projection of angular velocity on Z-axis 

The variables p, p, and r characterize the lateral motion of the 
airplane . 

By assuming the nondependence of the lateral motion of the air- 
plane on its longitudinal motion, the equations of the lateral dyna- 
mic stability can be derived. 

For steady rectilinear flight without sideslip 
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I = Y= Z *L=M = lSf = 0 
p = q = r = 0 

W SB V «i 0 

U = V 


( 1 ) 


where V la the Telocity of the airplane. 

After a certain lateral disturbance , the magnitudes p, r, and t 
will vary with time, whereas the magnitudes q, v, and u may be con- 
sidered constant. 

Because at the initial instant p=r=v=0, the disturbances 
of the magnitudes p, r, and v are denoted by the same symbols: 

dp = P 
dr = r 

dv = T 

Evidently, the lateral motion of the airplane will be characterized 
by three equations: (a) the equation of the equilibrium of the for- 

ces along the Y-axis, (b) the equation of the equilibrium of the 
moments of the forces about the X-axis, and (c) the equation of the 
equilibrium of the forces about the Z-axis. By the fundamental 
equation of mechanics, the product of the mass of the airplane by 
the projection of the absolute acceleration on the Y-axis is equal 
to the projection of all external forces acting on the airplane on 
the Y-axis . 

In the general case, the projection of the absolute accelera- 
tion of the Y-axis is equal to (reference 2) 

dv 

dt + ru - 

and the equation of equilibrium of the forces along the Y-axiB will 
he of the form 



I 


) 
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For the case under consideration! 

u «= V and w = 0 

In equation (2), Y denotes the sum of the projections on the Y-axis 
of all the external forces acting on the airplane so that 


Y 


+ *w 


(3) 


where 

Y a projection of aerodynamic force on Y-axis 

Yy projection of weight on Y-axis 

In what follows, it is assumed that the aerodynamic force acting 
on the airplane during a disturbed motion and the small deviations 
from the condition of steady rectilinear flight depend linearly on 
the disturbances, and because at the initial Instant no lateral force 
acted on the airplane the aerodynamic force Is given by 


v dY SY SY „ 


If the small magnitudes ^ p and ^ r are neglected, the fol- 
lowing equation is obtained: 


dY 


(4) 


In symmetrical flight there are no gravitational forces aoting on 
the airplane along the Y-axis. If the airplane rolls, a component 
appears on the Y-axis. When small oscillations occur, the force 
acting along the Y-axis for a rolling angle $ Is equal to 
(figs. 2 and 3) 


G- oos Q • $ 


I In a rigorous derivation, the term , pw drops out as a mag 
nitude of second order of smallness. 


\ 


I 
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■where 0 is the angle of inclination, of the flight path. Similarly, 
for the rotation of the airplane about the Z-axis hy an angle t 
the component along the Y-axis (fig. 4) is 

G sin 9t]/ 


Hence, for small deviations 

Y w = G cos 0$ + G sin 0*^ = G cos 0(<J> + tan 0tjr) 

But 


G cos 0 m Qx, ^L. S 


Therefore, 


2 

Y t = C L £L. S($ + tan 8*tf) 

C* 


(5) 


For small oscillations of the airplane, it may be assumed that 


v = - Y sin p a - Yp 

When this relation is taken into account and equations ( 3 ), (4), 
and (5) are substituted in equation (2) 

mY (~ ft + r ) " If $ + °L '^IT + tan ( 6 ) 

The equations of equilibrium of the moments about the X- andZ-axes 
are of the form (reference 3) 


A 




C 


dr 

dt 



N 


( 7 ) 


The moments of the aerodynamic forces L and IT on the free motion 2 
of the airplane after a disturbance will depend on the angular veloc- 
ities p and r and the angle of sideslip (3, For the small dis- 
turbances p, r, and p,. the aerodynamic moments may be represented 
as follows by a linear dependence on p, r, and p. 


2 The free motion of the airplane after a disturbance with locked 

controls , 
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L= ^ g+ ^ P+ ^? r 
w SN _ dfl ciN 

N ^ p + 35 p + 3? r 


Satiations (7) now "become 

.dp _ dr SL a &L dL 

A It 32 dt = 3F p + 35 p + 3? 


dr 

dt 


dp ■ SN 
dt 


SN 


Sn 




( 8 ) 


If equation (6) is differentiated together with equations (8), a 
system of three linear homogeneous differential equations Is 
obtained because p *= #/dt and r » d^/dt. 


mV 


[ c»> 

c dr dp ag |N dN 

C dt ' 35 at *3p p + p + I? r j 

Nondimenslonal magnitudes are now introduced by the equations 

* - <V *¥ s * 


L . Oj ^ Sb 


N 


2 

=nf 81 


(10) 


A a i 


mV 
"A 4 


* 

Gy is here the lateral force coefficient which should not be 
confused with the previously used notation for the lift coefficient. 
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SL 

3£ P 




dL pb - pY 2 

a J* 2T 
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r, dV^ 


p £ ££_ St, „ P l_ HZ- 

Sib 2 **P 2 

a 27 


Sb 


rl> 


dN - _. 

35 P - Pn p ^ 


dU 

3? r 


ra r 


Sc T 

7 b '# 

Scz 
Z P “ 3jT 


n P 


9c n 

w 

5cz 

s ££ 

27 


^3, 

5c, 

l *‘ h A 


6 27 


Sc n 

”p- 3E£ 


° 27 


Sc n 
= S £b 


° 27 


If all terms on the left in equation (9) are 
the nond imena iona 1 magnitudes sure introduced, the 


transferred and 
system becomes 
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_ d2p 

mV — £ - mV 
dt d 


27 dr 
b 4t + 


pY2 dp py2 27 _ 

~ 2 ~ Syp dt + C L ~ S» — (p + r tan 6) 


= 0 


*A 


mb 3 2V dg mb£ 27 dr pV^ 

4 b it" 1 ® 4 b dt' 2 




SbZpP 


0 


mb 2 2V dr mb 2 27 dp pT 2 _ pV 2 _ pY 2 
~~dt- i E'T"~dt-2 sb V “ 2 Shn r r ~ “ S *V 


0 


If the first equation of the system Is divided by p^V^S^/m, 

2 2 

the second by i A — T , and the third by i c °— 3 St> } the following 
equations are obtained: 

( m \ 2 d 2 p 1 dp m 2m _ 2m m dr 

pSY j dt 2 + 2 yp dt pSY + 2 pSb P ’ pSb pSY dt + 

C L 2m 

T 7§b tan 0Xr = 0 


m dp *p _ % m dr *r _ 

" i A 3 + pSY dt - i A P ' i A pSY dt -•* i A r * 0 


_ % _m_ d£ _m_ dr fr _ 

“i c P"ic P "ic pSV_dt + pSY dt " % r ■ 0 


By setting 


T = 


m 

pSV 

2m 

pSY 


( 11 ) 


and by introducing a nondimens ional time by the equation _ 


t = T 


( 12 ) 
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>5 

D 

-4 


Moreover, if the signs in the second and third equations are changed, 
the following system of equations in nondimsns ional form is obtained: 


d 2 g 

dt 2 


1 

+ 2 r P 





tan Sr = o) 


is 


p 


dp l v - ^ dr 

--f+v^P+T = 

dt *A. *A dt 



0 


> (13) 


*c 


P 



p 


% dp 

+ 1 

k dt 


dr 

dt 


°r - 
+ % r 


0 


J 


SCy 3c n 5c z 

The magnitudes y^ *= n^ *= and a are determined 

from the usual tests in the wind tunnel and are denoted as the sta- 
tic derivatives. The magnitudes I p , I r , n^, and n r characterize 

the moments of the aerodynamic forces in rotation and for this rea- 
son they are known as the rotary derivatives. 


2 . THEORETICAL DETERMINATION OF DERIVATIVES OF 
LATERAL DYNAMIC STABILITY 

As has teen pointed out, all aerodynamic coefficients are 
referred to a system of axes: the X-axis Is parallel to the veloc- 
ity and passes through the center of gravity of the airplane at the 
initial instant; the Y-axis passes along the span of the wing to 
the right of the pilot; and the Z-axis Is perpendicular to the first 
two and directed downward 3 , The positive direction of the angles of 
rotation and of the moments are assumed to he the following: 

about the X-axis, motion from the Y-axis to the Z-axis 

about the Y-axis, motion from the Z-axis to the X-axis 

about the Z-axis, motion from the X-axis to the Y-axis 

This system of axes is fixed In the airplane so that the X-axis 

coincides with the velocity of the airplane only at the Initial 


3 This system of axes is in accordance with the terminology 
adopted by CAHI, May 1939. 
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instant of time. For a disturbed motion of the airplane, the X-axis 
of the fixed system of axes makes an angle 6, which is termed the 
sideslip angle, with the flight velocity. 

The moments about the X- and Z-axis are denoted by L and If, 
respectively, and the angular velocities about _ the X- and Z-axis UJy 
p and r, respectively. 


ROTARY DERIVATIVE 2 p 

In the following discussion, all the derivatives in the belcw- 
stall range will be considered. In the case of rotation of the air- 
plane about the X-axis (angular velocity p), an aerodynamic damping 
moment arises about the same axis. This moment may be determined by 
the equation 


L . 


_ 2k 


dV 2 


2Y *P 2 


Sb 


where 

Zp moment coefficient 

^ = p nondimens ional velocity of roll 

Here as usual (reference 1), a linear dependence of the aerodynamic 
reaction on the disturbance is assumed. The rolling moment is pro- 
portional to the fourth power of the linear dimension, and therefore 
the moment on the tail surface for an angular velocity p is neg- 
ligibly small as compared with the moment acting on the wing. 

Hence, in computing Z p only the moment acting on the wing shall 

be taken into consideration. 

On rotation of the airplane about the longitudinal X-axis with 
angular velocity p > 0 (the right wing drops); at a flight veloc- . 
ity Y, it is found that the angle of attack increases on the right 
wing. At the distance y from the plane of symmetry 

A a » tan a = 
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Similarly , on the left side of the wing the angle of attack decreases. 
A moment therefore appears opposing the rotation. Thus, for p > 0, 

L < Oj and the coefficient of the derivative of the rolling moment 
■with respect to the angular velocity of roll is always (at be low - 
stall angles of attack) negative (fig. 5) Zp < 0. Figure 6 gives 

curves by which Zp for wings of various aspect ratios and tapers 
may be determined (references 4 and 5). 

Example: r\ = 1.94, X «= 6.34. 

From figure 6: 1? = -0.465. 

It can be seen from the curves that Zp increases in absolute 
value with an increasing aspect ratio of the wings and that Zp 
decreases with increasing taper Zp 

The value of Zp is computed on the basis of the general theory 

of the wing with account taken of the change in the distribution of 
the circulation along the span due to rotation without considering 
the change in velocity along the span. At below-stalling angles, 
the derivative Z p theoretically does not depend on the initial 

twist of the wing nor on the angle of attack of the entire wing. 

The magnitude Zp varies little for different wings, as may be seen 
In figure 6. The mean value is = -0.45. The sign of Zp may 

change only in the above- stalling range In the region of autorotation 
of the wing. 


ROTARY DERIVATIVE np 

In the rotation of the airplane about the X-axis, the angle of 
attack varies along the wing span and as a result the lift force 
and the induced drag will not be the same at symmetrical elements. 
Hence, when the airplane rotates about the X-axis, the airplane will 
be acted on by a moment of the aerodynamic forces both about the 
X-axis (Zp) and about the Z-axis (np). As before, consider that 

the moment about the Z-axis arising for angular velocity p is pro- 
portional to the angular velocity of rotation, so that 

If = H n p qSb 

where q. is the dynamic pressure. 
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The moment of yaw due to the angular velocity defends little on 
the tail surface so that in the computations' for the dynamic stabil- 
ity only the moment due to the wing is taken into account. 

For positive rotation in roll p > 0 (right side of wing drops), 
the results show that on elements of the right and left wing the lift 
force vector L and the vector of induced drag D are rotated and 
changed (fig. 7). The airplane will therefore be acted on by a 
moment about the Z-axis (moment of yaw), which for negative angles 
may be positive and with increase in the angle of attack of the wing 
varies linearly so that at moderate angles of attack np is always 
negative (np < 0 (fig. 7)). Curves of the values np/c ^ (refer- 
ence 4) for untwisted wings 4 of various tapers and aspect ratios are 
shown in figure 8. The curves give the values of rip /a^ where 

is the angle of attack in degrees and is computed from the line of 
zero lift. Hence, in order to actually compute the value of np, it 

is necessary to multiply the value obtained from the curves for the 
given wing (given aspect ratio and taper) by the angle of attack (in 
degrees) computed from the line of zero lift. 

For example, .q = 1.94, A = 6.34, e 9° (c^ is computed from 

the angle of zero lift). From figure 8, rip/c^ * -0.00395 and, there- 
fore, (np)a = 9° = -0.00395X9 = -0.0356. 

ROTARY DERIVATIVE I r 

For positive angular velocity r, the left wing moves forward, 
that is, the left wing has a greater velocity than the right wing. 
Therefore, the lift force on an element of the left wing will be 
greater than the lift force on the symmetrical element of the right 
wing and as a result the airplane will be acted on by an aerodynamic 
moment tending to lift the left wing, that Is, a positive moment 
about the X-axis. With increase in the angular velocity r, the 
moment L Increases in magnitude and the derivative of this moment 
with respect to the angular velocity r = rb/2V will be positive 
(at positive C^) 5 . 


4 The value of np little affects the characteristic of the 
dynamic stability and therefore a correction for the twist of the 
wing need not be introduced. 

Here, as previously, linear dependence of the moment on the 
angular velocity Is assumed as is permissible because only small 
angular velocities are considered and the terms containing the 
squares of these velocities may be neglected. 
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The coefficient Z r is 
tional to the angle of attack eta. Evidently , "both the wing taper 
and the aspect ratio will affect the value of I r . Curves of the 
variation of for untwisted wings of various tapers and 

aspect ratios are given in figure 10. From the curves it may he 
seen that l r decreases with increase in the taper and increases 

with increase in the aspect ratio. 

It is also apparent in figure 10 that * r /a& for various types 

of airplane changes little with wings of different tapers and aspect 
ratios. The mean value is \ T fa^ = 0.017 

0.018 

In order to obtain the value of Z r at a given angle of attack, 
the value Zp/a^ obtained from the curves must be multiplied by the 

angle of attaok (in degrees) for which the computation is conducted 
(the angle of attack is measured from the line of zero lift) . 

For example, t] = 1.94, A = 6.34, = 9°. From figure 10 

Zr/a^ = 0.0185j Z r = 0.0185X9 = 0.1665. The vertical tail surface 
will also affect the derivative Z r of the airplane. The part of 
l r due to the tail may be obtained from the equation 

*r = a t (if) sin 2( x - a*) , (14) 

where 

a t = (dCj/dcOt (fig. 11). 

Kj. empirical interference coefficient at tail, mean 

value = 0.8 

Z, distance from center of gravity of airplane to 

z rudder hinge 


rb 

— gY" 2 r qSb 
l r > 0 

determined for the wing and is propor- 
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Stf area of vertical tall surface with part of fuselage 

b wing span 

S wing area 

x angle between line of zero lift and line connecting 

center of gravity of airplane with geometric cen- 
ter of vertical tail surface (including part of 
fuselage) (for a low wing airplane x > 0) 

<% angle of attack (from line of zero lift), degrees 


Example : 04. 

x 


2.05 

12 ° 


Kfc . 0.8 


h 

b 


0.389 



0.093 


To find 2 r ^ ’ for =9°, from equation 14 T 


I rt = 2.05 X0. 8(0. 389) a X 0.093 X 0.105 = 0.00242 


CORRECTION ON DERIVATIVE l r FOR TWIST OF WING 

In order to Introduce a correction for twist, it is necessary 
to use figure 12. Only the practical .determination of the correc- 
tion will be considered here®. The correction for twist for wings 
of various aspect ratios and the taper T| = 3 can be found by use 
of figure 12. 

For a taper different from 3, the correction would be somewhat 
different but the error will be very small and it may therefore be 
assumed that the correction for twist does not depend on the wing 
taper. Assume, for example, the correction for the wing is found 
to be A = 6.34, the twist of the wing starts at the distance 
0.226 b/2 from the root, increases linearly, and at the tip attains 
the value 1.5°. 

At the distance y = 0.226 b/2 from the root, the twist Is 
equal to zero (fig. 13(b)). In figure 13(a) at y = 0.226 b/2, 

A-l = 0.001. Subtracting this value from D (fig. 13(a)), 0.0168 is 


6 This method and the curves are taken from reference 5, where 
a more detailed description of the curves is found. 
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obtained. This number is plotted in figure 13(c) and corresponds to 
an angle of twist equal to zero; therefore , on the curve of fig- 
ure 13(c) the point (0, 0.0168) is plotted. It can be seen that at 
the distance 0.4 b/2 the angle of twist is 0.3° (fig. 13(b)) and 
from figure 13(a) Ag = 0-003; 2 r /a = D - Ag = 0.0158. In fig- 
ure 13(c) 0.3° is plotted as the abscissa and -0.0158 as the ordi- 
nate, and so forth. A certain area is then obtained on the plani- 
metering for which the correction in Z r is obtained. 

For the given case, l v = 0.0144. 

This correction does not depend on the angle of attack of the 
wing 7 . 


For the total value of 2„ 

r 


*r 



“a + *rt + ^r twist 


(15) 


Example: find 2 r : 7\ = 6.34, tj = 1.94, *= 9°, a^. = 2.05, 

K t = 0.8, ij b = 0.389, S tf /S = 0.093, * = 12° 


The twist starts at 0.226 of the half span, increases linearly, and 
at the tip reaches the value 1.5° 

2 r - {l r /cL & )a & + 2 ri . + 2 r = 0.1665 + 0.00242 + 0.0144 . 0.1833 


ROTARY DERIVATIVE n^ 

During an angular velocity of rotation about the Z-axis (angu- 
lar velocity of yaw r) a yawing moment arises, which depends on 
the vertical tail surface, the fuselage, and the wing. At an angu- 
lar velocity r, on account of the difference in the velocities at 
the symmetrical elements of the wing, different drags are obtained 
as a result of which there arises a moment about the Z-axis (moment 
of yaw N) . Evidently, the moment N due to the wing will be neg- 
ative for positive values of r (fig. 14). Similarly, for the 
angular velocity r > 0 a negative yawing moment is found to arise 
because of the fuselage and the tail. 


7 The value of 2 r ^. y is positive for the case where the 

angle of attack at the tip of the wing is greater than the angle of 
attack at the root section, . 
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The yawing moment will increase in absolute value with increase 
in r; hence, 


N = 2y Sh 
Hj. < 0 (fig. 14). 

The effect of the tail and the wing on the derivative n r will now 
he considered. 


Effect of Vertical Tall Surface on n^. 

In the oscillation of the airplane about the Z-axis, that is, at 
an angular velocity r, the tail surface is found to be acted upon by 
the force 



rlt 

V 2 


cos( x - a ) 
8 > 


where x is’ the angle between the line of zero lift and the line 
connecting the center of gravity with the geometric center of the 
vertical tail surface with part of the fuselage. For a low-wing 
airplane, the angle x > o. In the expression for the force during 
oscillations about the Z-axis there enters cos (x-Og.) because the 

system of axes Is fixed. If the derivative of the moment of this 
foroe is taken about the center of gravity with respect to the non- 
dimensional angular velocity r = rb/2V, 


Hi 


rt 


-2 ' 


/dC L \ pS tf v2 t 


\da } t b 


cos^f x - a^) 


Finally, the nondimens lonal derivative rip Is obtained by dividing 
this expression by pSV^b/2 


■»- «?)>&)' 

In this equation no account Is taken of the effect of the fuselage 
on the flow at the vertical tail surface. The indeterminancy of 
the magnitude of the interference at the vertical tail, surface makes 
the theoretical computation of the moments due to the vertical tail 
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surface inaccurate. .An empirical coefficient Kj. characterizing 

the velocity interference at the tail surface is therefore intro- 
duced. Moreover } setting 


gives 


"rt 



- a t 



(* - eta) 


By permitting a small error } it may he assumed that cop^(x 
and then 


n rt 




°a) s#1 


(16) 


where 


distance from center of gravity of airplane to rudder hinge 
h wing Bpan 

S-^f area of vertical tail surface with part of fuselage 
S area of wing 

For the coefficient K^, the value 0.8 may he taken. 

The value of a-fc depends on the aspect ratio of the vertical 
tail surface. 


Figure 11 shows the dependence of a t on the experimentally 
obtained aspect ratio (reference 6). The mean value of the 

coefficient a^ -for airplanes of the usual type is 

= 2.2 

In order to compute the aspect -ratio of the vertical tail surface 
Aj., the span of the tail b^ must he taken as shown In figure 11. 

In the tail surface area there is included the end of the 
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fuselage (hatched area in fig. ll). The effect of the fuselage is 
not taken into account in equation (16). In the same manner as for 
determining the damping of longitudinal oscillations, the coeffi- 
cient 1.25 in equation (16) is introduced to take account of the 
fuselage . Then 


Example : 


n rtf 


2 . 5K+a-jj 




K-f. — 0.8 as 2 . 05 


2 t /b sr 0.389 


S-^f/S as 0.093 


(16a) 


By formula (16a), 

n rtf *» 2. 5x 0.8x2. 05X (0.389) 2 X 0.093 = -0.0577 


Effect of Wing on Derivative n^. 

For positive angular velocity about the Z-axis (r > 0) , the 
velocity on the left wing increases and on the right decreases. 
Hence, the drag on the left wing will be greater than on the right 
wing and there appears a moment about the Z-axis acting against 
the rotation (fig. 14). 

The damping moment of from the wing Is due to the induced 

and profile drags. 

Curves for the determination of nj, of the wing due to the 
induced drag are given in figure 15. In order to determine n^, 
from figure 15, it is necessary to multiply the value of n r /a^ 2 
by Oq 2 where is the angle of attack of the wing computed 

from the line of zero lift in degrees. The part n^. of -the wing, 

which is obtained as the result of the profile drag,' is small and 
need not be taken into account®. 


®The curves for computing this part of the derivative^ n r are 
given in reference. 5. 
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Example. Find, of the wing: T| = 1.94, A = 6.34, = 9°. 

In figure 15, n r /a a 2 = -0.000147, n^ = -0.000147X 9 2 = -0.0119. 

The total value of the derivative n^. for the airplane may he 
obtained from the formula 

«r ■ ^tf + n rw ( 16b ) 

Example: Find n^, of the airplane = 9^. 

The characteristics of the wing are: T) = 1.94 arxL A = 6.34. 

The characteristics of the tail surface are Stf/S = 0.093, 

7 t /h = 0.389, = 1.14. 

n^. *= n rtf + = -0.0577 - 0.0119 = -0.0696 

It must be borne in mind that at times the interference of the 
fuselage has a strong effect on n^ . For example, for one model a 

value of rij. of the tail and the fuselage was obtained equal to a 
third of the derivative n^, for the tail only. Hence, a sufficiently 
reliable determination of n r can be obtained only by experiment. 

The theoretical equations for determining n^, may in individual cases 

give considerable error and should be considered only as rough 
approxlmat ions ® . 

Equations have been derived herein for the computation of the 
rotary derivatives, that is, the coefficients of the moments obtained 
during rotational motion of the airplane (yaw or roll). The deter- 
mination of the coefficients of the forces and the moment of the 
aerodynamic forces In sideslip follows. These coefficients are called 
static derivatives. 


%arly In 1940 at CAHI, work will be completed on the determina- 
tion of the derivative for 10 to 15 models. From this investi- 

gation, the value of the derivative n r can be obtained for' similar 
airplanes. 
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DETERMINATION OP STATIC DERIVATIVE Z p 

For a positive sideslip angle 6 the left wing moves forward, 
the flow under the left wing tending to lift it, that is, for 
P > 0 the rolling moment L is positive and 

L = LpP t= pZpcjSb > 0 
Zp > 0 (fig. 16) 

Experiment shows the static derivative Ip may he found by 
the formula 

l p alrp^ *= 0 ^ + * + l p fus + l pt (1?) 

X = 0 ^ 

where 

iff angle of dihedral of wing 

X sweepback angle of the wing 

The characteristics of the individual terms in the equation 
will now be discussed. The magnitude gives Ip for the 

X*0 

wing in the absence of sweepback and dihedral. The derivative 
(2p)^ _ 0 arises from the change in distribution of the circulation 

X =0 

over the wing span as the result o;f the sideslip. Test curves show 
that this derivative increases with increasing angle of attack. 

This magnitude cannot be computed on account of the absence 
of test data. The plan form of the wing considerably affects the 
value of this coefficient. A rectangular wing without rounded tips 
gives a value much greater than a wing with rounded tips (reference 7). 

The effeot of the dihedral of the wing is taken care of by the 
Slo 

term f 3?ke an 6l° of wing dihedral V' i 0 nieasured between the 

plane of the chords and the plane perpendicular to the plane of 
symmetry of the airplane and passing through the root chord (fig. . 17 ) • 
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Experiment allows that 2p is proportional to the dihedral ^ 
with a coefficient of proportionality cilp/cty, which may he deter- 
mined from figure 18, which gives the values of the derivate Ip 
for a wing without a straight midwing section. 

Examplej Find lg of the wing: p = 1,94, 7v= 6.34, ifr = 4.25° 

According to figure 18, 5lp/SV f = 0.0133. 

Slg 

Ip -g— «= 4.25X0.0133 = 0.0565 

For a wing with a midwing section it is necessary to subtract from 
the value obtained from figure 18 the amount 

(^)mw - O-O 2 (t^ 

where b^/b is the ratio of the span of the 

span of the wing. The angle f in this caBe 
shown in figure 17(b). 

Example: T| >= 1.94, Tv a 6^54, y = 4.25°, 

Zp = 4.25 [o/0133 - 0.02(0.226) 2 ] 

Tests show that the tip of the wing greatly affects the value 
of the derivative 2p. The wing tip in vertical projection may be 

such as shown in figure 19. 

From specially conducted tests (reference 7 ) , it was found 
that if the wing tip is of the shape shown in figure 19(c), the 
effective dihedral is increased by 1 to 1.5°. 


(18) 

midwing section to the 
is to be measured as 

b^/b = 0.226 
= 0 . 0523 . 


Effect of Sweepback 

The angle of sweepback is measured between the lines passing 
through the fooal line at 0;25 chord from the leading edge of the 
wing and the plane perpendicular to the axis of the fuselage. The 
sweepback angle X will be considered positive for the direction 
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shown in figure 20. For a positive sweepback at a positive sideslip 
angle , the flow about the left wing is improved and the lift force 
on the left wing is increased. As a result, an increase in the roll- 
ing moment is obtained for positive sweepback - " 


X >0, 2 p > 0 

The effect of sweepback is taken into account in equation (17) by 

' 9*6 . — 

the term X At the present time sufficient data, by use of 

which the effect of the sweepback for different wings oould be 
determined, are unavailable. F. Veick (reference 8) recommends 
for rectilinear wings of aspect ratio 6 the equation 



0.0045 C L 


(19) 


where X is the sweepback angle in degrees. 

Example. X <= 2.5°, = 0.73(a = 9°) 

From equation (19) 

i 

= 2.5(0.0045X0.73) = 0.00822 

Effect of Position of Wing Relative to Fuselage 

The previously described relations give Zp for isolated wings. 

Tests show, however, a very great effect of the position of the wing 
with respect to the fuselage on the derivative Zp. It was found, 

for example, that a center wing gives a value for Zp approximately 

agreeing with the value for the isolated wing. A high wing increases 
the effective dihedral by 1° to 1.5° and a low wing decreases the 
effective dihedral by 2° to 5°. This change in the effective dihed- 
ral occurs, however, at small and medium angles of attack. At large 
angles of attack different results may be obtained, which at the 
present time, due to the absence of sufficient test data, cannot be 
even approximately predicted. 


On account of the resistance of the flow, all projections at 
the lower surface of the wing increase Zp and conversely all pro- 
jections on the upper surface decrease 2p. In a fixed system of 
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axes, the X-axis passes below the center of pressure of the vertical 
tail surface j therefore, the Tertical tail surface -will also give a 
rolling moment in sideslip. If the angle between the line connecting 
the center of gravity with the center of pressure of the vertical 
tail surface and the velocity is denoted by - x 

I{3t - ^t a t s~ a ^ n ( x " °a) (20) 

Example. ELj. = O.S, l-j^/b *s 0.389, S^^/S e 0.093, a^. s 2.05, x 1 * 12®, 

% - 9° 

By equation 20 

Zp t = 0.8X0.389X0.093X2.05X0.0523 = 0.00597 


It follows from the foregoing discussion that the derivative 

l a may be determined with sufficient accuracy only by experiment. 

P 

These equations for determining the effect of the dihedral and 
sweepback of the wing may be used for estimating the changes In the 
airplane parameters after tunnel tests . For example , If from com- 
putation of the test data It is found that it is necessary to 
increase Zp by 0.026 and the parameters of the wing are q = 3, 

\ «= 7 (without mid wing) the required increase in the dihedral 
is obtained from the relation tsjf X 0.0133 = 0.026 (the value 
0.0133 is obtained from fig. 18) whence A^ * 2°. If In another 
case it is necessary (for example, from considerations of longi- 
tudinal stability) to decrease the sweepback by X = 5®, the 
change in Zg can be found. (The computation is conducted for 

the Initial state of flight of the given airplane). Let = 0.6 



0.0045C L = 0.0045X0.6 


0.0027 


Alp = 0.0027X5 = 0.0135 


Hence, decreasing the dihedral by 5°, Zg decreases by AZg = 0.0135. 
In order to compensate for the decrease in Zg the dihedral of the 
wing must be increased 


AtjrXO . 0133 = 0.0135 
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whence 


= 1 ° 

that is, the dihedral muBt he increased by 1°. 


DETERMINATION OF STATIC DERIVATIVE 

The static derivative np = SCn/Sp characterizes the change 
in the moment of yaw with change in the sideslip angle and depends 
mainly on the area and the shape of the vertical tail surface and 
the fuselage. The center of pressure of the fuselage for the usual 
arrangement is located ahead of the center of gravity and for this 
reason the yawing moment due to the fuselage for p > 0 is posi- 
tive, that is, the fuselage is unstable. As may be seen from fig- 
ure 21, np < 0. 

The magnitude of the coefficient np is determined analyti- 
cally by the equation 


where 


np = np f + np t = K 


ff If h ftf 

S b “ E t a t b S 


( 21 ) 


lateral area of fuselage 
If length of fuselage 

The coefficient Kp may be found from figure 22 (taken from 
reference 9) 

where 

Zf length of fuselage (positive) 

Znf distance of center of gravity from nose of airplane (positive) 

h maximum height of fuselage (positive) 

The coefficient Kp(Kp > 0) may be determined if the ratios 
Z^/h and are known. 
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Example: l^/h = 

From figure 22 it can be seen that Kp = 0.106. 

If S^/S = 0.35 and 2f/b = 0.595 then 
npf = 0.106 X 0.35 X 0.595 = 0.0221. Equation (26) ie suit- 
able for airplanes with unstable fuselages. Sometimes the fuse- 
lage is stable (as is rarely the case) and this eq.uation cannot 
then be applied. The value of a^. may be obtained from figure 11. 

The mean value for is 0.8. 

Example: 2 t /b = 0.389, S^/S = 0.093, \j. = 1.14. 

From figure 11, a^ = 2.05 

' n pt = “ 0-8X0-389X0.093X2.05 = - 0.0593 

np s= Upp + Up-{; = 0.0221 — 0.0593 = — 0.0372 

Let the area of the tail surface be increased such that 
AS tf /S e 0.04. Then for the increase in the coefficient np 

Anp = - 0.8X0.389X0.04X2.05 = - 0.0268 

Thus, when it is required to increase the vertical tail area, the 
span must be increased but not the chord. If the vert leal tail 
surface is Increased by increasing the chord, decreases and 

therefore a^ also decreases; as a result for such an increase in 
S tf the weathercock stability (np) is found to undergo only 
slight increases or does not change at all. 

The equations for determining the value of np must be consi- 
dered to be only approximate. For a more accurate determination of 
np, recourse must be made to experiment. As in the case of deter- 
mining the derivative 2g , these approximate equations may be used 
for computing the derivative np after the tail area is increased 
or decreased. 
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Effect of Propeller 

In the present Invest igation, the effect of the propeller slip- 
stream on the, dynamic -stability derivatives Is not considered. The 
derivatives n^ and n^, evidently will be greatly affected by the 

propeller slipstream. Up to "the present time, however, sufficiently 
accurate test data, by which a correction for the effect of the pro- 
peller on the tail surface may be Introduced, are unavailable. 


DETERMINATION OF STATIC DERIVATIVE y p 

The static derivative yg (fig. 23) is always positive 
(yp > 0) . The derivative yg characterizes the increase In the 

lateral force (y) acting on the airplane during an increase in 
the sideslip angle (p). Evidently, the value of yg depends 

essentially on the size and shape of the fuselage and of the ver- 
tical tail surface. Diehl (reference 1 0) recommends the following 
empirical equation for computing the derivative yg: 



b2 f 





y p « 0.12 -g- 


(22) 

- 

Example: 

h = 14.53, = S.B7, - S = 33.4 



t 


y 3 = 0.12 blf - 0.12 14-53 X 9.67 = 0#453 
p S 33.4 



__ 


For airplanes with usual fuselages, the mean value of yg Is 0.4; 
therefore, the derivative yg has little effect cn the character- 
istics of dynamic stability. 


3. EXPERIMENTAL DETERMINATION OF DYNAMIC-STABILITY 

DERIVATIVES 

The dynamic-stability derivatives are divided into two groups; 
namely, the derivatives of the moments and forces with respect tc 
the angular velocity, which are called the rotary derivatives, and 
the derivatives of the moments and forces with respect to the side- 
slip angle, which are the static derivatives. 
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The rotary derivatives may he determined in a -wind tunnel on 
special apparatus. The procedure of these tests is, however, rather 
complicated and they will not be considered herein (reference 11). 

The static derivatives n^, Zp, and yp may be readily 

obtained from the results of the usual tests in the wind tunnel. 

For this purpose it is necessary to have the coefficients of 
yaw and roll and the lateral force (C n , Cj, and C y ) for vari- 

ous angles p at the given angle of attack. Curves are given In 
figures 24 uo 26 for the change in the coefficients C n , Cj, and C y 

with sideslip angle. 

0 

If these coefficients are given in the fixed system of axes 
assumed at CAHI (fig. l), then in determining the derivative no, 
the tangent to the curve C n = f (6) for 6=0 must be drawn 

and the slope of the angle of inclination measured (taking account 
of the scale). Special attention must be paid to the determination 
of the slope of the curve C n = f(p). The slope of this curve 

characterizes the weathercock stability. For stable airplanes, the 
curve is like that shown in figure 24, that is, the coefficient C n 

decreases with increase in the sideslip angle p. For neutral or 
laterally unstable airplanes, the curves differ to a marked degree 
from thoee in figure 24 near 6 = 0 only, in the range from 6 = -3° 
to p s +3°, where for neutral airplanes these curves run parallel 
to the axis of abscissas and for unstable airplanes the slope of the. 
curve changes sign, that is, the coefficient C n in this range 

increases with increase in p. For a sideslip angle less than 
B = -3° or greater than B = +3°, the curves in all cases are as 
shown in figures 24, that is, the coefficient C n decreases with . 

increase in the sideslip angle p. Hence, in determining the 
weathercock stability, the points lying within the range p = +3° 
must be considered. During the tests, points must be taken every 
l c angle near p = 0 but beyond jp| =3° points may be taken 
every 2° to 3°. The computations required for obtaining the deriv- 
ative are presented in figures 24 to 26. The coefficient 57.3 is 
introduced into the computation for converting degrees into radians 
in measuring the angle p. For a sufficiently accurate determina- 
tion of the derivatives np, Zp, and yp, it is necessary to com- 
pute them for 4° or 5° angles of attack. (See appendix 1.) 


4 
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Determination of Moments of Inertia of Airplane 

In order to compute the lateral dynamic stability of the air- 
plane, the coefficients of the moments of inertia, which can be 
obtained from the following formulas, must be known: 

4A 4C 

1a — r • o and ip = p 
A mb ' 2 0 mb 41 


where 

A polar moment of inertia with respect to X-axis 

C polar moment of inertia with respect to Z-axis 

b wing span 
m mass of airplane 

Because for various conditions of flight of the airplane (dif- 
ferent angles of attack a) the axes X and Z rotate with respect 

to the airplane, - it is evident that the moments of inertia will vary 
likewise. If the moments of inertia of the airplane about the prin- 
cipal axes (A’ and C) are known, the moments of inertia about the 
X- and Z-axes may be computed by the following equations: 

A =A ' cos 4 * <P + C ’ sin 2 cp y C = C * cos 2 9 + A 1 sin 4 * cp, E « — — sin 2<p 


(23) 

Computations show that up to angles of attack a = 15° to 20° it 
may be assumed with a sufficient degree of accuracy that A. s A • 
and C = C’, that is, the moments' of inertia about the X- and 
Z-axes are equal to the corresponding moments of inertia about the 
principal axes. As will be shown, the centrifugal moment of Inertia 
has little effect on the stability characteristics and therefore E 
need not be computed. 

At the present time no sufficiently accurate statistical for- 
mulas exist by which the moments of inertia of an airplane may be 
quickly determined and it is therefore necessary to resort to 
analytical computation In each case. Formulas and the computing 
procedure of the moments of Inertia are given in reference 12, 
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The values of 
C.049 < i A < 0.089 


i^ and 1 q change -within the ranges 
and 0.09 < i q < 0.154. 


Determination of "Airplane Density" p 
The "airplane density" 

d 

where 

G/S load/sq. meter, (kg/m 2 

y unit weight of air 

b span, meters 

The computation of the dynamic stability is usually conducted for 
the altitude H = 3000 meters and therefore y = 0.908 kilograms 
per cubic meter. For the altitude H =3000 meters for standard 
atmosphere 

G 

p = 2.20-f- (24) 

Example : 


p is determined by formula (ll) 

2 - 
2m S 

= oSb ~ 7 b 


— = 78 kilograms per square meter 
s 


b = 14.53 meters 


p = 2 .2 


78 

14.53 


11.8 





V- 


is a ncnd imens ional magnitude 
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Determination of t 

The magnitude t is determined by equation (11) 

t _ 2 L_ 
dSV 

In determining m for H = 3000 meters 


1 1 - 2 - 


(25) 


where G/S is in kilograms per square meter , Y is in meters per 
second, and t is in seconds. 

Example: G/S = 78 kilograms per square meter and V =47.4 meters 

per second 


t = 1.1 = 1.1 ~~ = 1.79 sec. 


t- 

o 

CVJ 


4. SOLUTION OF SYSTEM OF EQUATIONS OF MOTION 

A solution of the system of equations (13) will now be sought 
In the form 

(3 = Ee^; p = Pe** and r = He^ 

By substituting In equations (13) and dividing by e^, the following 
equations are obtained: 


(AS + -|ygA^B + p np P + p tan 9 

(v) B ’ ( x '^) p + (^ X4 ^) R -° 

% b + (5 4 ^ x ) p -( x -^) r '° 


R = 0 



MCA TM 1264 


35 


The system of linear homogeneous equations has a solution if the 
determinant of the system is equal to zero. 



A 





( 


Ct 

-g- tan 6 



b 

±A 


-A 





*0 




In expanding the determinant, an .equation of the fourth degree 
in A is obtained, which is called the characteristic equation. 

A^A" + AjA^ + Ag\^ + A^A + Aq = 0 (26) 


where 
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Analysis of the coefficients shows that up to angles of attack from 
15° to 20° the magnitude 1^, has little effect on the roots of the 

characteristic equation and without large error i^ = 0. 

If in addition, this analysis is restricted to the considera- 
tion of level flight, that is, tan 0=0, the following equation is 
obtained: 


A 4 = 1 



0 ‘ 4 2 V l A ’■O ' ^ ) 


Because 

y 3 > 0, J s > 0, tip < 0 
Zp < 0, Z r > 0, np < 0, n r < 0, C L > 0 
n > 0, i A > 0, i c > 0 
The following relation is usually true: 

a 3 > °> a 2 > °> A 1 > 0 

The coefficient A Q , the free term of the characteristic, equation, 
may be either positive or negative. 


2074 - 
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Example : 


jyp.O.Zi; ^ = 1.175; 


= - 0.312 


h 


7.27; = 3.11 


• 7 ^- = - 0.52; 
In 


cl 

-gr = 0.37; 


°r 

T“ — ~ 0.76 


4 = 11.8 


In order to find the coefficients of the characteristic equation, 
use is made of equations (27): 

A 3 " I + IjJ “ °- 24 ' [ ' 7 - 27 * 0.7e] " S' 27 

2 \*a h ic *a J V 2 y B/ L*a *cJ 1 g 

* {7.27X0.76 + 0.52X3.11} - 0.24 {- 8.03] + 11.8X0.312 = 12.7526 


"P 


°L *£ 


i + ^ 9 i 
X C w A 


/l \ *p "r Z r1 ^ ha l & 

= 2 y P ' i I iT ' f + U T ~ i " 

V 7 |^A 1 C 1 C ^-Aj x c x a x a 

> 0.24 {7.145} + 11.3X0.312X7.27 + 11.8X1.175X0.52 + 

11.8X0.37X1.175 = 40.809 

C-r /!_ rig 2 q n^A 

A 0 - 1* T ^ j - ^ (- 3.1W0.312 + 


1.175X0.76) • - 0.3362 
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If the roots of equation (26) are denoted by A^ Ag, ^3, and 
the solution of the system (13) is obtained in. the form 

A-i t Aot , A/t 

6 = B^e + Sge + B36 + B^e 


Aj-t A^t A^t A 4 t 

p = Pj© + Pg® + -^3® + ^4 e 


A^t Agt A jt A^t 

r cr R^e 4* RgS ^ P3© ^4® 

Evidently, in order that the sideslip angle and the angular 
velocities of the airplane 8, p, and r should decrease with time 
(that is, in order that the airplane should be dynamically stable), 
it is necessary that all real roots be negative and that the com- 
plex roots have negative real parts. In the general case, the solu- 
tion of an algebraic equation of the fourth degree is Laborious. In 
the given case, however, by making use of the special characteristics 
of the equation of lateral stability, this equation may be quickly 
and 3 imply solved. The determination of the .roots of equation (26) 
takes no more than 20 minutes. 

For airplanes of the usual type, the characteristic equation 
of the lateral stability has two real roots of which one is very 
large and the other very small and two complex conjugate, roots. The 
small real root is denoted by A^, the large real root by Ag, and 

the complex conjugate roots by A3 and A 4 . Because A^ is very 
small, the following equation may be written with a large degree of 
accuracy; 

A^A+Aq = 0 

from vhioh A-^ is immediately determined^ 0 

•^Practical, computation shows that the value A = Mq/a^ is 
accurate. For example, if A ^ is graphically determined, the 
obtained value always agrees with A^^ = -Aq/a^. 
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* - J 2 (28) 

Now It can "be seen that If Aq < 0, then A-^ > 0 and the airplane 
is unstable . 

Example.: A^ = 40.809 Aq « -0.3362 

From formula (28) 


*1 



-0.3562 

40.809 


0.00825 


In determining Ag, must he recalled that A g is large and there- 
fore if the powers of Ag below the third are neglected, the approxi- 
mate value is obtained 


^ 2 4 + A 3^2 3 w 0 

whence 


~ - A 3 


Computations on many airplanes show that the best approximation for 
Ag is the value^ 



(29) 


Because i A > 0 and at below- stalling angles of attack < 0, 

Ag < 0. 

By setting 

A (A) * k ^ +■ AjAS + Ag\2 + A^ A+ Aq 

I 


II Usually the value of Ip/ i^ gives a very small error in the 
determination of the root Ag. However, a more accurate determination 

Is necessary to compute the complex roots of the characteristic equa- 
tion (more accurately, the real part). 
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for a more accurate determination of the root A is plotted 

as a function of A near a value equal to A ~2p/i A , by laying 

off A on the abscissa and A on the ordinate. The intersection 
of the curve A(\) vith the abscissa gives the accurate value of 
A o . Because 


A* (A) e= 4A 3 + 3A3A 2 + 2A2^ + A]_ 

and as a polynomial for large numbers has the sign of its greatest 
term and Ag is negative, usually A'(A 2 ) < 0. Hence, if after 

the first trial the point falls above the X-axis, the next value 
of A g in smaller absolute value is taken, whereas if the point 

falls below the X-axis the next value in greater absolute value is 
taken. Usually, three points are sufficient for drawing the curve 
A (A) and determining A 2 . 

Example: Find Ag from the equation 

A 4 + 8.27 A 3 + 12.75 A 2 + 40.809 A - 0.3562 = 0 

Because l^/i^ * -7.26, the value Ag 1 «‘-7.2 is used for the first 

value of Ag » , 

The order of computation is clear from table I. 

The second, third, and fourth powers of A, which are required 
for the computation, may be found in appendix 2, which gives the 
powers of numbers from 2 to 15 in steps of 0.2 and covers the entire 
range of values of A encountered in computation. The curve A (A) 
from which Ag = -7.3 is taken, is given in figure 27. After having 

found A^ and Ag, the conjugate complex roots A^, A 4 are obtained 
from the following considerations: 

Equation (26) may be put in the form 
A^ + AjA 3 + AgA 2 + AjA + Aq = ( A - A^)(A - Ag)$^ + aA + b) * 0 


or 


A 4 + [a -(Aj^ +Ag2A 3 + [b - a(A^ + Ag) + A^Ag] + 
[aAjAg — b(A^ + Ag)] A+ bA^Ag = 0 


2074 



By equating the coefficients of the same powers of A the following 
equation is obtained: 


Aj r a. •( A^ + Ag ) 

Ag =' h - a(A^ + ^g) + 

A 1 = aA L A 2 " + ^g) 

A q = bXjXg 

whence b = Aq/A^\,j but by (28) A^ = -A^/A^, hence 

^i 

1=-^ (30) 

2 

The value of a is found from. 

Ag = b — a (A -^ + A g ) + A^Ag 


If a and b 
the equation 


a 


b + A-jAg ■“ Ag 

A, + A, 


'1 - r, '2 

are known, the roots A 3 and A^ 


(31) 

are obtained from 


A^ + aA + b * 0 

^3,4 = - f ± 1 ^ 

If As = i + it) , A 4 = I - irj and 12 


if 


-*- 2 These formulas are written on the assumption that b > a 2 /4; 
h < a 2 /4, the formulas become 



and in this case the airplane is unstable. If A-^ 0, then b ^ 0 
and instability results. A change in the sign of A^ is generally 

due to the Iosb of weathercock stability of the airplane. 
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t «= 


a 


2 





(32) 

(33) 


Example : Find the complex roots of the equation. 

A. 4 + 8.27\ S + 12.75 A? + 40.809 A- 0.3362 . 0 
It has "been established that 

A-^ * + 0.00825 and A2 = - 7,3 

An png 

From equation (30), b = - = 5.61 and by equations (31) 

b +^2 - A 2 5>61 + o.00825(-7.3) -12.75 „ _ 

a * -"A-TT^ — o:dd825 -'7.3 " °* 990 

It can be seen with the aid of equations (32) and (33) that 


l - -ft 


T] « 


0.990 


0.495 


>\J b - (|) 2 = 2.315 


5. MOTION CHARACTERIZED BY ROOTS K lt A g , A 3 , AND A 4 

The disturbed motion of an airplane is character iied by the 
expressions 


0 = b/ 1 * + Bj.®* + B 3 a^ + B 4 , Mt 


„ - ^ _ AS* _ A4 t 

p s P-^e + Pge + P^e + P^e 


A jt 


- _ \l* _ A2* _ A3* _ A4 t 

r = R^e + Rg© + Rj© + R^e 
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This motion may he broken down into three components 


A _ t t A - 1 

I p-^ = B^e j Pp = PpS J rp = Bp® 


A„t _ \gt _ At 

II p 2 = ®2 e 5 ^2 = P 2 e J r 2 = B 2 e 


A_t A,t _ A_t A,t _ A_t A.. 

Ill p 3 « B 3 e ° + B 4 e , p 3 . P 3 e 3 + P 4 e & , r 3 . Rge 5 + H 4 e * 

The computation of the coefficients B^, Bg, B 3 , B 4 , P^, Pg, and so 

forth will not he considered in detail. It should he noted that of 
all these magnitudes only four are arbitrary and are determined from 
the initial conditions, whereas the remainder can he obtained from 
these four magnitudes. 

Prom a numerical analysis of these coefficients the following 
conclusions may he made: 

I. Consider the group 


Ai t 

Pp = Bp® 


_ Ait 

r^ = R^e 

The magnitude P^ is small in comparison with B-j_ and R^. The 
root A^ is "very small. Evidently, if the airplane is stable, 

A^ < 0 and the magnitudes p (sideslip), p (angular velocity of 

roll), and r (angular velocity of yaw) decrease with time (with 
increase in ^ = t/r ) . 
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The root characterizes spiral stability of the airplane. 

In order to clarify the physical significance of spiral stability, 
the motion after a roll of the airplane is considered. Let the 
airplane roll in such a manner that the right wing drops. Whereas 
in the horizontal position of the wing the lift force balanced the 
force of gravity, the force of gravity will now give a projection 
on the Y-axis and the airplane begins to sideslip (p) on the 
right (lowered) wing. As a result of the sideslip, two aerodynamic 
moments appear: the moment ZpXp tending to rotate the wing in 

the horizontal position and the moment npXP decreasing the side- 
slip angle 3. Because of the action of the moment n^XP the air- 
plane begins to yaw, that Is, an angular velocity r appears, which 
gives rise to two new moments: the moment n^r opposing the moment 

npXP and a moment Z r xr, which opposes the. moment ZpXp. 

Thus, the characteristics 2g and n^ act favorably in elim- 
inating the disturbance and the characteristics n^ and l r act 
unfavorably (with respect to spiral stability). 

As can be seen from equation (28) spiral stability Is charac- 
terized by A q, which is determined by the equation 

°L 
f* T 

A 0 ~ i A Iq ^r 1 ^ “ ®r l p) 

where 

C L 

H > 0; “g- > 0 

l r > 0; np < 0; n^. < Oj Zp > 0 

Hence, there will be stability for the case where r^Zp 
value is greater than Z r np whence 


z p 

> 

*r 

*0 


iV 


in absolute 
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The degree of stability is characterized by the time in which the 
initial disturbance decreases by one-half 


e^L* 



J^t = l n | = - 0.693 

but 

t-i 

^spir = " 0.693 ^ (34) 

In the case of instability .X-^ > 0 and the time of decrease of the 

initial disturbance by one-half will be negative. In the case of 
stability X^ < 0 and t ^ will be positive 

T =1.83 sec X = + 0.00825 

1.83 

^spir * * 0.693 0.00825 ” " 154 sec 

It has been found that 2.5 minutes after the initial instant the 
initial disturbance will increase to twice the value. Such motion 
evidently will not be observed by the pilot. 

Experience shows that spiral instability is not felt by the 
pilot if the time of increase of the initial disturbance to twice 
the value is large. It is to be supposed that if 

l^spir) > 40 - 50 sec 

at the Initial flight conditions the airplane will be satisfactory 
with regards to spiral stability. 

II. Consider the group 

Xpt 

8 2 = B 2 e 
— Xgt 

P2 = P 2® 


Xgt 

r 2 = K 2 e 
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Numerical analysis shows that during a disturbance Bg and Eg are 
small in comparison with Pg; that is, this group is found to char- 
acterize a motion which differs little from a pure roll. - It has 
been shown that 


and because at below- stalling angles 2^ < 0, p will decrease very 
quickly with time because \g is large. Generally, at the below- 
stalling angles this motion is not noticeable after 1 second. 

III. Consider the group ■=- 


= B 3 e 




+ B 4 e 


\ 4 t 






+ P.e 
4 


k 4 t 


Because and 

X 3 = i + in; 
A 4 « t ~ in; 


A,t 

B 5 e + B 4 



B-a 




X 4 t 

e 


\,.t A.t 

r 3=V +S 4 9 

B 3 and B 4 are conjugate complex quantities 

: bj + ib 4 as B(oos 3 + i sin 3) * Be 
= b^ - ib 4 = B(cos 3 - i sin 3 ) = Be 1 "^ 


= Be 


i3 (£+ in)t -ip {I - iT))t 

e + Be e 


= Be U [e (T!£ + P)i + e" (llt+P)i ] 
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But. e ( 1 l t + P) i + e" ( T l t+ P) i = 2 cos(T)t + fi) j and, therefore, if 2B = B 


Similarly, 


M 3 ^ 4 * -St 
= B.^e + B^e = Be cos(tit + p) 


_ A^t A.t _ £t _ 
p = .P_e + P.e = Pe cos(r|t+y) 


_ X._t \ 4 t _ £t _ 

r 3 = + ®4 e = cos(rit + B) 

This group characterizes the oscillatory damping of the motion. 

The time required- for decreasing the amplitude by one -half can 
now te found. 

£t 2 1 


£t 2 = - 0.693 


to = 


Therefore 


Example: T = 1.83 


to = ~ 0.693 r 

t 

£ = - 0.495 


tg = - 0.693 -r^ 


= 2.56 sec 


The period of oscillation is determined from the relations 

tjT = 2n 

_ T 
T - t 


T = 2rt 1 

n 
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Example: T = 1.83 sec T) = 2.31 

T . 6.28 -|^|| = 4.98 sec 

With regard to the oscillatory stability of the airplane, the 
following condition must be imposed: The airplane should be stable 

in oscillation and the time of decrease of the initial disturbance 
by one-half (tg) should be less than 14-1.5 times the' oscillation 

period 

t g < 14-1.5T 

The stability in oscillation is affected mainly by the lateral static 
stability n^. The other coefficients also affect the oscillation of 

the airplane but to a smaller extent. Whereas an airplane may have a 
small static instability, it may nevertheless be dynamically stable 
but the period of oscillation (T) and the time of decrease of the 
initial disturbance by one-half (tg) will be large. At greater 

static instability, the airplane will be unstable in oscillation; 
that is, the amplitude of oscillation of such an airplane will 
increase with time. For a coefficient of lateral static stability 
of the order = 0.01 and greater, the least disturbance will 

terd to increase the roll of the airplane (a kind of slow spin). 
Whereas from tests in the wind tunnel of the model with locked con- 
trols it was found that the airplane is neutral or even has a small 
weathercock stability, in flight such an airplane will be unstable 
because of the mobility of the control surface. For this reason, 
from the requirement of oscillatory stability it follows that the 
airplane should be statically stable over the entire range of flight 
angles of attack and especially at below-stalling angles of attack. 


6. ORDER OF COMPUTATION OF DYNAMIC STABILITY 
Computation of Airplane Northrop 2E for Dynamic Stability 

In this section, the entire dynamic-stability computation is 
presented from setting up tables of the structural parameters to 
obtaining the final dynamic stability characteristics. At the 
start of the computation, a table must be set up in which all para- 
meters of the airplane required for computing the lateral dynamic 
stability are included. 
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I. G-, v eight of airplane, kilograms 

II. Wing 

1. b, wing span, meters 

2. b^, midwing- sect ion span, meters 

3. S, wing area, square meters 

4. Cq, chord at root, meters 

5. c , chord at tip, meters 

"C 

6. V, dihedral angle - angle "between plane of chords and. 

plane perpendicular to plane of symmetry of airplane 
and passing through the chord at the tip, degrees 
(fig. 15) 

7. X, sweptback angle - angle between lines passing through 

focal line taken at one-quarter chord from leading edge 
of wing and plane perpendicular to axis of fuselage 
(X > 0 for forward sweep! ack) , degrees (fig. 19) 

8. 7, angle of twist of wing, degrees 

III. Vertical tail surface 

1. s tf ^ ver ' fc I° a l "bail surface area (area of tail with part 

of fuselage, (fig. 11), square meters 

2. b^., tail span (measurement shown in fig. 11), meters 

3. I distance from center of gravity to hinge of rudder, 

meters 

4. Angle between line of zero lift and line joining center 

of gravity of airplane with geometric center of verti- 
cal tail surface (with part of fuselage) 

IV. Fuselage 

1. , area of projection of fuselage on plane of symmetry of 

airplane, square meters 



so 
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2. Z f , length of fuselage, meters 

3. h, maximum height of projection of fuselage on plane of 

symmetry, meters 

4. Z^f, distance of center of . gravity from nose of fuselage^ 

meters 

The magnitudes computed under I, II, III, and IV" are always 
positive. 

The table of structural parameters for the Northrop airplane is: 


G = 

2600 kilograms 

S tf 

=3,11 square meters 


S = 

33.4 square meters 

b t 

= 1.88 meters 


b = 

14.53 meters 

h 

= 5.65 meters 


^mv = 

3.29 meters 

Sf 

= 11.69 square meters 


c 0 " 

2.9 meters 

*f 

= 8.67 meters 


°K - 

1.5 meters 

h 

= 1.45 meter b 

w 

If = 

7° 30' 

2 nf 

= 2.11 meters 

n 

X = 

2° 30’ 

X 

= 12° 



7 = 1° 30* 

The airplane is a low-wing single- engine type 


|-78 kg/m 2 


S. 


tf 


= 0.093 


n - c^ = 1,94 


= 0.389 


_£ 

s 


0.35 


v £2 
K = S 


6.34 ^ = 1.14 


= 0.226 




G 


= 2.2 s = 11.8 



0.244 
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Determination of Aerodynamic Coefficients 

1. The curve of variation of the coefficient of lift (C^) 

against the angle of attack (a) obtained from tests in the wind 
tunnel for a model of the Northrop airplane is given in figure 2. 

The angle of attack is measured from the wing chord. Because C^ 

must he determined for angles of attack measured from the line of 
zero-lift, table II is obtained. 

2. The rotary derivatives Zp, Z r , n^, and n^. were determined 

from the aforementioned equations and curves and also experimentally 
on a special apparatus (reference 11), 

(a) Determination of Zp 

The derivative Zp is obtained from the data of the table of 

structural parameters (fig. 6). The curves experimentally obtained 
and by figure 6 are given in figure 30, Theoretically, the deriva- 
tive Zp does not depend on the angle of attack but during a test 

a certain change in Z with change in a & is usually obtained. 

The agreement obtained, as shown in figure 30, is sufficiently good 
and for computing the dynamic stability the value of Zp may be 

taken from figure 6 without resorting to experiment. 

• (b) Determination of n^ 

The theoretical value of Up is found from figure 8 by using 

the given t] and k in the table of structural parameters. The 
value of np/ag. taken from figure 8 must be multiplied by 

(angle of attack In deg. from the zero lift line). The values of 
n^ obtained from test and from figure 8 are given in figure 31. 

It is obvious that a rather large disagreement results. A syste- 
matic numerical analysis of the dynamic stability equations shows, 
however, that the dynamic stability characteristics depend little 
on the value of np and for this reason only an approximate value 

of Up need be known. Hence, for computing the dynamic stability, 

the value, which is obtained from figure 8, may be used with suffi- 
cient accuracy. 

(c) Determination of Z r 
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The determination of the derivative 2 r must be divided into 
the following three stages: 

(1) Determination of 2 r for an untwisted wing 

(2) Determination of the correction on 2 r due to the twist 

of the wing 

(3) Determination of the part of Z r due to the vertical tail 

surface 

1. The derivative Z r is obtained from figure 10 by the given 
T) and A. of the wing. The value 2 r / taken from figure 10 must 
be multiplied by a. (from the line of zero lift). The curve of 

2 r against for the Northrop wing- is presented in figure 32. 

The dashed curve gives Z r for the wing without twist. 

2. As the Northrop wing is aerodynamically twisted such that 
the angle of attack at the tip of the wing is greater than the 
angle of attack at the tip of the midwing by 1.5° > figure 12 must 
be used to determine the correction. This correction evidently 
does not depend on the angle of attack and therefore the entire 
straight line undergoes parallel displacement. The correction for 
twist is 


J p, twist = 0*0144 

3. The correction on Z r due to the vertical tail- surface is 
obtained from equation (14). The required value of a-j. is deter- 
mined from. figure 11. The continuous line in figure 32 shows the 
derivative 2 r of the Northrop airplane with account taken of the 

twist of the wing and of the vertical tail surface. The figure 
also gives the points obtained, from experiment. It is apparent that 
the agreement is good and that the derivative Z r may be computed 

theoretically without recourse to experiment. 

(d) Determination of 

The derivative n^. is made up of a part due to the fuselage 

with vertical tail surface and a part due to the wing. The part 
n^ due to the fuselage with tail is computed from equation (16a). 


I VLOZ 
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The value a^ is obtained from figure 11. The part n^. due to the 

■wing is computed from the data presented in figure 15. In order to 
obtain r^, the value r^/a^ from figure 15 must be multiplied by 

the square of the angle of attack. The angle of attack must be in 
degrees and measured from the zero-lift line. The sum of these two 
parts gives the derivative xi^ of the airplane. The values of the 

derivative n^. for the Northrop airplane obtained from test and. 
from computation are shown. in figure 33. 

It should be remembered that the derivative n~. plays a great 

part in computing lateral dynamic stability and, for this reason it 
is desirable that n^, be computed as accurately as possible. 

Remembering also that equation (16a) may in certain cases give a 
large error, this magnitude must be experimentally determined. 

(e) Static derivatives yp, Z p, np 

The static derivatives are determined from standard wind-tunnel 
tests; therefore these magnitudes as a rule must be obtained as shown 
in section 3. The values of yp obtained by equation (22) (curve) 

and by experiment (points) are presented in figure 34. Equation (22) 
generally gives a sufficiently accurate value of yp. The values 

of Zp for the Northrop airplane obtained from test and from compu- 
tation are given in figure 35. The upper dashed curve gives Zp 

for .the value V = 7.25° (dihedral angle). In order to take account 
geometrically of the effect of the fuselage for a low wing in the 
formula, the dihedral must be decreased, as has been pointed out, by 
2° to 5°. An airplane is considered the chassis of which is not 
streamlined and all projecting parts on the lower surface increase 
the effective dihedral; therefore, for the computation, the ^ angle 
minus 3° is used, that is, an angle of 4,25°. The lower dashed curve 
in figure 35 gives the value of Zp obtained from computation 

- 4.25°). The value dZp/dty is obtained for the given p and 
A from figure 18. 

The effect of the midwing on the derivative is taken into 
account by the decrease in dz^/Sty by an amount that nay be obtained 

by use of equation (18). The center line in figure 35 gives the 
value of Zp with account taken of the part due to the vertical tail 

surface and the continuous line gives the value Zp, if the sweep- 
back of the wing is considered. 
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’ The values of* *r 5c n /d|3 obtained from test and by equa- 
tion (21) are given in figure 36. In equation (21) , the value 
Kj. = 0.8 is used; the value of the coefficient is obtained 

from figure 22 and from figure 11. 

Notwithstanding the sufficiently good agreement that was 
obtained for the Northrop airplane by computation and test for the 
values Ip = and n^ = c)C n /d{J, these formulas must be con- 

sidered only approximate and In all cases the values Ip and np 
determined from the test curves C n = f(6) and Cj = f ^ ( p ) . 

Moments of Inertia of the Northrop Airplane 

The moments of inertia, as has been previously stated, may be 
obtained by computation. For the Northrop airplane, the moments 
of inertia are determined by full-scale experiments according to 
the method explained In the report by U. A. Pobyedonostsev (ref- 
erence 13). The results of the investigation are as follows: 

I A = 0.0579, i G = 0.0959 

where the values are given in nond imens ional form. 

For convenience of computation, a table of all coefficients 
that enter the characteristic equation must be set up as shown In 
table 3, in which for the derivatives the test values were taken 
and the angle was computed from the line of zero-lift. 

From the data of table 3, the coefficients of the character- 
istic equation A^, Ag, A 1 , and A Q are computed by equations (27). 

The values of these coefficients for the Northrop airplane computed 
for the angles of attack = 1°, 5°, 9°, and 13° are given in 

table 4 and the values A-^, Ag, £, and r|, computed by the equa- 
tions of section 4, are given _in figures 37 and 38. A comparison 
of the computed results with flight tests of reference 14 Is made 
in figure 39. 

Complete agreement was obtained for the period of oscillation T. 
The disagreement that exists for tg (the time of decrease of the 

amplitude to l/2) may be explained by the fact that no account wag 
taken in the computation of the effect of the propeller slipstream 


20^4 
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on the vertical tail surface. The magnitudes A^ and Ag in 

flight have not teen determined; thus it is impossible to check 
the accuracy of the computations of these values. It can be seen 
from figure 39 that the computation gives sufficiently good results. 


7 . APPROXIMATE FORMULAS 

In this section, equations are given for estimating roots 


These formulas 
X l> X 2> and 


quickly from the given aerodynamic coefficients. 

may also be used in computing the corrections on a^, Ag, S, and q 

if after computing the lateral dynamic stability it is decided to 
change the area of the vertical tail surface or the dihedral of 
the wing. 

For the root Ag, the equation Ag = ^p/*A can ** e use( ^» The 
approximate expression for Ap, which determines the spiral stabil- 
ity, is obtained from A-^ = -Aq/A-^ and because 


A 1 


X A 


2 in ) + U U 


W *C 


therefore 


Ai = 


(37) 


2 ^r n p) 

ic - + Vp 

The time of decrease of the initial disturbance to l/2 is obtained 
by the equation 


^spir = ' °- 693 7^ 

For the magnitude i (giving the time of decrease In the ampli- 
tude of oscillation tg to l/2) the following equation can be ' 

used: 


l = 


fr 

^ - 2 


A n 


2 


(38) 
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It can "be seen from equation (38) that the value cf 4 increases 
in absolute value with increase in magnitude of the damping of the 
yaw Up ' and the value y^, which characterizes the lateral force 

in sideslip. The magnitude of 4 drops in absolute value with 
increase in 1q of the nondimens ional radius of inertia about the 

Z-axis. Vibrational instability will occur for 4 > 0 and, as oan 
be seen from equation (38), vibrational instability may be obtained 
for a large spiral stability if 0, which is usually obtained 

for large values of the derivative 2p and in the absence of 
weathercock stability, that Is, when >0. For this reason, 

when the airplane has small weathercock stability or instability, 
the wings should not be given a large transverse dihedral because 
the derivative 2^ will then be large and this condition will lead 

to.excessive spiral stability and, what is particularly important, ■. 
may give rise to oscillatory instability. 

If the area of the vertical tail surface is changed, the change 
in rip is obtained by the equation 


ft t\ 2 AS tf 

An r ~ " 1,oa t ytTJ s ' 

where 

a^ = (dC^/ScOf. (obtained from fig. 11) 

2-t distance from hinge of rudder to center of gravity 

AS-fc change In vertical tail surface area (AS > 0 for increase 
In area; AS < 0 for a decrease) 

The change in 4 is obtained from the equation A4»An r /2iC, 

The new value of .. 4 is obtained from 4 * =4 + A 4 . From the newly 
obtained value 4 ' , tg ie from equation (35). 

%2 - “ 0»693 7*7 ‘ 

where T = 1.1 (G/S)/V (for an altitude of 3000 m) . 
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Because A^_ may Be approximately determined by the equation 



and 


I p A x 

^ i n ‘ w ^ lere = " ~ 


the following equation is obtained: 


/ r b /°L M ^3 
~ [“ ^p \ 2 " ic J " 


(39) 


whence the period T is determined by the equation 


T = 6.28 £ where t 
T| 


= 1.1 


G 

S_ 

V 


The greater r) is the smaller the period of oscillation is 
found to be and therefore the period decreases with decrease in the 
absolute value of the derivative np, with increase in Iq, the non- 

dimensional radius of inertia with respect to the Z-axis, and with 
decrease in the value of the derivative Zp. If the expression under 

the square root sign becomes negative (as may occur in the case where 
np > 0, that is, In the absence of weathercock stability) the magni- 
tude q does not characterize the period of oscillation but together 
with the magnitude % determines the asymptotic deviation of the 
airplane from Its course. 

The previously derived approximate equations give very good 
agreement with exact computation for the Northrop 2-E airplane. In 
other cases, generally Bpeaking, the agreement nay be less complete 
but a qualitatively correct result will be obtained^. 


l^At the time the present article was being prepared for pub- 
lication the previously derived formulas were checked for a large 
number of airplanes. In almost all cases the approximate formulas 
give very good agreement with the exact formulas and they may there- 
fore be used in computing the lateral dynamic stability of the 
a irplane . 
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APPENDIX 1 

Up to the present time Investigations in the -wind tunnel have 
been conducted In a system of axes different from that assumed in 
this paper and for this reason equations are presented herein with 
the aid of which the old. coefficients my he reduced to the coef- 
ficients used in this paper. 

The sideslip angle 3 (previously y) in the diagrams of the 
wind-tunnel tests has the same sign as 3 assumed in the report. 
The moment coefficients have the opposite sign. Moreover, up to 
this time the moment coefficients were referred to pV 2 in order 
to reduce them to nondime ns ional coefficients. If the yawing 
moment (in the old notation C^) Is based on the distanoe between 

the center' of gravity and the rudder hinge and the rolling moment 
(C^) is based on the span. 


Cl 





Cy c — 2Cjr 


If the results from the tunnel investigations (C m , C^, C z ) 
are given in a gravity system of coordinates , 

and Sc z /S3 are obtained from the curves and the coefficients Ip, 
np, and yp are obtained from the equations 

SCjjjx 

H c = " ~W X 57 * 3 


5c_ 1+ 3C~. 3c z 

n 3 -3T-- 2 b"^r* 57 - 3 ; ^3 - 2 W* 57 - 5 


If the results are given in a system of chord axes, all moments 
and forces are based on oY 2 , the rolling moment C’^ is based, on 

the span, and the moment of . yaw C’^ Is based on l t (the distance 
between the center of gravity and the hinge of the rudder). 


C =0’ 
mx mx 


sin a; C = C 1 
my my my 


*t 

cos a - ~ C ’ 


mx 


C. = C* 
z z 


^t 

cos a + ~ C’ 


sin a 


7.CHA. 
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j 

j 


and finally 

3C 2 

z p ■ W 


^ l ^mx 

= - 2 ~w 57,3 = 


Sc 

~3p 


mx 2 t ^my 

cos a + 5p“ sin a 


^n 

n P “ '3[T 


2 

' _c 


SC' 


t ^ 'my mx . . 

to" ~w~ 003 a ’ ~5p sin a 


57.3 


Sc, 


Sc 


= w = ■ 2 “sr x57 - 3 


57.3 
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SECOND, THIRD AND FOURTH POWERS OF NUMBERS FROM 2.0 TO 15.0 ° 


A 

A 2 


A 4 

A 

A 2 

A 3 

A 4 

A 

A 2 

A 3 

A 4 

2.0 

4.00 

8.000 

16.000 

6.6 

43.56 

287.5 

1897.5 

11.0 

121.0 

1331.0 

14641.0 

2.2 

4.84 

10.65 

23.43 

6.8 

46.24 

314.4 

2138.1 

11.2 

125.4 

1405.0 

15735.2 

2.4 

5.76 

13.82 

33.18 





11.4 

130.0 

1481.5 

16890.0 

2.6 

6.76 

17.58 

45.70 

7.0 

49.00 

343.0 

2401 

11.6 

134.5 

1561.0 

18106.4 

2.9 

7.84 

21.95 

61.47 

7.2 

51.84 

373.2 

2687.4 

11.8 

139.2 

1643.0 

19387.8 





7.4 

54.76 

405.2 

2998.7 





3.0 

9.00 

27.00 

81.00 

7.6 

57,76 

439,0 

3336.2 

12.0 

144.0 

1728.0 

20736.0 

3.2 

10.24 

32.77 

104.36 

7.8 

60.84 

474.5 

3701.5 

12.2 

148.8 

1816.0 

22153.3 

3.4 

11.56 

39.30 

133.63 





12.4 

153.8 

1906.6 

23642.1 

3.6 

12.96 

46.66 

168.0 

8.0 

64.00 

512.0 

4096 

12.6 

158.8 

2000.4 

25117.4 

3.8 

14.44 

54.87 

208.5 

8.2 

67.24 

551.4 

4521.2 

12.8 

163.8 

2097.0 

26843.5 





874 

70.56 

59277 

4978.7 





4.0 

16.00 

64.00 

256.0 

8.6 

73,96 

636.0 

5470.1 

13.0 

169.0 

2197.0 

28561.0 

4.2 

17.64 

74.09 

311.2 

8.8 

77.44 

681.5 

5996.9 

13.2 

174.2 

2300.0 

30350.6 

4.4 

19.36 

85.18 

374.8 





13.4 

179.6 

2406.1 

32256.0 

4.6 

21.16 

97.34 

447.7 

9.0 

81.00 

729.0 

6561 

13.6 

185.0 

2515.5 

34210.2 

4.8 

23.04 

110.6 

530.8 

9.2 

84.64 

778.7 

7163.9 

13.8 

190.4 

2628.0 

36267.4 





9.4 

88.36 

830.6 

7807.5 





5.0 

25.000 

125.0 

625.0 

9.6 

92.16 

884.7 

8493.5 

14.0 

196.0 

2744.0 

38416,0 

r™ n 

o,c 

in it ni 

140.6 

rrrrl rt 

/OX. 2 

9.8 

96.04 

941.2 

922377 

1472 

■2UI76 

206373 

40658,7 

5.4 

29.16 

157.5 

850.3 





14.4 

207,4 

2986.0 

42998.2 

5.6 

31.36 

175.6 

983.4 

10.0 

100.0 

1000.0 

10000 

14.6 

213.2 

3112.1 

45437.2 

5.8 

33.64 

195.1 

1131.6 

10.2 

104.0 

1061.2 

10824.3 

14.8 

219.0 

3242.0 

47978.6 





10.4 

108.1 

1125.0 

11690.6 





6.0 

36.00 

216.0 

1296 

10.6 

112.4 

1191,0 

12624.8 

15.0 

225.0 

3375.0 

50625.0 

6.2 

38.44 

238.3 

1477.6 

10.8 

116.6 

1259,7 

13604.9 





6.4 

40.96 

262.1 

1677.9 
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TABLE I 




SBEB 



n 

X'” = 

-7.3 

n 

mm 

X n 

An* 1 

x ,m 

.tin 

A n A 

I 

A n X ,,,n 

4 

i 

2687.38 

2687.38 

2998.68 ' 

2998.66 

2839.8 

2839.8 

3 

8.27 

-373.248 

-3086.76 

-405.224 

-3351.2 

-389.02 

-3217.2 

2 

12.75 

51.84 

660.96 

54.76 

698.19 

53.29 

679.4 

1 

40.809 

-7.2 

-293.82 

-7.4 

-301.98 

-7.3 

-297.9 

0 

-0.3362 

1 

-0.34 

1 

-0.34 

1 

-0.34 




-32.58 


43.33 


3.76 


TABLE II 


a °a 

i° 

50 

90 

130 

B 

0.092 

0.44 

0.74 

1.04 


zone. 
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STABLE HI - COEFFICIENTS FOR COMPUTING 


DYMMIC STABILITY OF NORTHROP AIRPLANE 



°a = 1° 

= 5° 

0& - 9° 

= 13° 

yp 

0.53 

0.53 

0.48 

0.45 

l p 

0.068 

0.068 

0.068 

0.068 


-0.033 

-0.037 

-0.03 

-0.032 


-0.45 

-0.45 

-0.42 

-0.38 

i r 

0.042 

0.115 

0.180 

0.26 

“p 

-0.019 

-0.034 

-0.05 

-0.065 - 

% 

-0.055 

-0.0675 

-0.073 

-0.069 


1 


1 



% 

=17.3; 

i=- = 10.4 
H3 



TABLE Hla 



o 

i — i 

n 

«P 

il 

oi 

o 

% « 9° 

o 

to 

i — ( 

it 

c L 

2 

0.046 

0.22 

0.37 

0.52 

Is 

2 

0.265 

0.265 

0.24 

0.225 

Z P 

Ta 

1.175 

1.175 

1.175 

1.175 

*C 

-0.396 

-0.385 

-0.312 

-0.333 

h 

A 

-7.3 

1 

• 

CD 

-7.27 

-6.57 

h 





l Y 





*A 

' 0.727 

1.99 

3.11 

4.5 

fE 

■0.197 

-0.354 

-.052 

-0.676 

n r 




-0.717 

*C 

-0.572 

» 

-0.702 

-0.76 



Ji = 11. a 
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TABLE IV 


\An 


Q 

n 

A 0 

1 

8.637 

11.5 

41.04 

0.209 

5 

8.747 

12.89 

44.43 

0.1524 

9 

8.27 

12.75 

40.809 

-0.3362 

13 

7.51 

13.32 

44.14 

-4.025 
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Figure 20.- - Sweepbaok of wing. 



Figure 21. - DeriTatiye Bp. 


Figure 22. - Curves for determining effect of fuselage on derivative n. 



Figure 23. - Derivative y ft . 




£U/4 


NACA TM 1264 





76 


MCA 'EM 1264 



Figure 26. - Variation of lateral force coefficient C y vith ■ sideslip angle (3. 

dc y 0.0825*57.3 _ 

— io — “°- 5 



Figure 27. - Determination of root X 2 of characteristic equation. 
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Figure 29, - Curves of lift force and drag against angle of attack. 
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Figure 30. - Derivative 2 p of Northrop airplane. 
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0 5 a 0 10 15 

a 

Figure 32. - Derivative 7 r of Northrop airplane. 



Figure 33. - Derivative of Northrop airplane. 
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Figure 34. - Derivative y 
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Without decrease In \|/ 

\J; decreased "by 3° (effect of fuselage) 

Zp with effect- of tail 


Zp of airplane (computed) 




Figure 35. - Derivative Z 
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• • • Experimental j 

-- ■■■ ■ Computed from equation (21) 


Figure 36. - Derivative n^ = Sc n /Sp . 








Figure 37. - Values of roots of characteristic equation a 1 emd a 2 
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